Statistical Mechanics of Quantum-Classical Systems with Holonomic Constraints 
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The statistical mechanics of quantum-classical systems with holonomic constraints is formulated 
rigorously by unifying the classical Dirac bracket and the quantum-classical bracket in matrix form. 
The resulting Dirac quantum-classical theory, which conserves the holonomic constraints exactly, 
is then used to formulate time evolution and statistical mechanics. The correct momentum-jump 
approximation for constrained system arises naturally from this formalism. Finally, in analogy 
with what was found in the classical case, it is shown that the rigorous linear response function of 
constrained quantum-classical systems contains non-trivial additional terms which are absent in the 
response of unconstrained systems. 



I. INTRODUCTION 

Quantum-classical formalisms U are computationally 
useful approximations of full quantum mechanics. In the 
past few years, a mathematical formalism, which permits 
to derive surface-hopping schemes from full quantum me- 
chanics by means of a number of controlled approxima- 
tions, has been proposed 0- Most notably a consistent 
statistical theory of quantum-classical systems has been 
first introduced and then used in order to formulate 
the theory of nonadiabatic rate constants 0. Numer- 
ical calculations of quantum-classical rate constants in 
nonadiabatic chemical reactions have been performed for 
some model systems which are relevant for condensed 
phase Despite some numerical problems with long- 

time integration of surface-hopping trajectories, this ap- 
proach holds promise for the study of more realistic sys- 
tems with few quantum degrees of freedom but many 
classical particles. As a matter of fact, just recently, 
Hanna and Kapral ^ reported the results of a quantum- 
classical calculation of the nonadiabatic rate constant for 
the proton transfer process occurring in a two-atom com- 
plex immersed in a classical bath of diatomic molecules, 
which were represented by means of cartesian coordinates 
and holonomic constraints. This example shows how, 
within the quantum-classical formalism of Ref. 0, one 
can easily model classical baths which are, in principle, 
as complex as the state of art atomistic representation 
of a protein by means of holonomic constraints and force 
fields (for an example see Refs. [^). However, in order 
to do so rigorously, one must generalize the quantum- 
classical formalism of Ref. so that it can describe 
classical baths with holonomic constraints. This issue 
was not addressed in Ref. 6]. Thus, the motivation of 
the present paper is to formulate a consistent statistical 
mechanics of quantum-classical systems when holonomic 
constraints are employed to model the classical bath. 
Linear response theory for quantum-classical constrained 
systems is also given because it is particularly relevant for 
the calculation of nonadiabatic rate constants. 
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A rigorous formulation of linear response theory for 
fully classical systems with holonomic constraints has 
been introduced only recently by means of a re- 
interpretation of the formalism originally developed by 
Dirac B. This formalism |^ , based on a generaliza- 
tion [lOllll| of the symple ctic structure of classical brack- 
ets in phase space 0,13,^3' showed that unusual terms 
may appear in the response function of a classical system 
with holonomic constraints. In light of this result, the ex- 
tension of the theory of Refs. ^ to the quantum-classical 
case is a subtle issue. The identification of the symplectic 
structure of quantum commutators TEj and the general- 
ization of this structure to introduce a non-Hamiltonian 
quantum bracket Il5||. of which the quantum-classical 
bracket of Refs. 0, Ql is a particular realization, has 
also been proposed just recently and used to introduce 
quantum-classical Nose- Hover dynamics . 

In the present work, the non-Hamiltonian commutator 
of Ref. Tsj and the classical Dirac bracket of Refs. |^ are 
combined in order to formulate the statistical theory of 
quantum-classical systems with holonomic constraints. A 
Dirac quantum-classical bracket, which exactly conserves 
constraints, is easily introduced and used consistently 
to define the dynamics and the statistical mechanics of 
quantum-classical systems with holonomic constraints. 
In particular, it will be shown how the momentum-jump 
approximation must be modified in order to consider cor- 
rectly the back-reaction of the quantum degrees of free- 
dom on the constrained classical momenta. The station- 
ary constrained quantum-classical density matrix will be 
derived and linear response theory formulated. As al- 
ready observed in the classical case 0, non-trivial terms, 
associated to the action of the perturbation on the phase 
space measure of the unperturbed constrained system 
and the phase space compressibility introduced by the 
perturbation itself, must be considered in general. It re- 
sults from the derivation that these additional terms are 
zero if quantum-classical variables, to which the pertur- 
bation is coupled, depend only from particlepositions. In 
this case, the rate formulas derived in Refs can be 

applied to a constrained systems with the only additional 
requirement of using the correct constrained stationary 
density matrix. However, different perturbations could 
require to evaluate all the terms in the response function 
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of quantum-classical constrained systems and one must 
be aware of their existence. 

If the work presented in this paper is considered to- 
gether with that of Ref. where non-Hamiltonian 
quantum commutators where introduced and used to in- 
troduce quantum-classical Nose dynamics, it is readily 
realized that a unified formalism for defining generalized 
dynamics and statistical mechanics in quantum-classical 
systems is now available. There are reasonable expec- 
tations that this formalism could be used in the future 
in order to attack the problem of long-time numerical 
integration of quantum-classical dynamics. 

This paper is organized as follows. In Sec Dirac 
formalism, as given in Ref. |^, is shortly summarized. 
In Sec. mil the results of Ref. [Ts'l are quickly pre- 
sented by showing the generalized symplectic structure 
of quantum-classical brackets, which are a particular 
realization of non-Hamiltonian quantum commutators. 
In Sec. IIVI the matrix structure of quantum-classical 
bracket, is used to combine it with Dirac bracket. In 
such a way the Dirac quantum-classical formalism is in- 
troduced and equations of motion which preserves the 
constraints exactly are given. In Sec. Ivl such equations of 
motion are represented in the adiabatic basis and the cor- 
rect momentum-jump approximation for a constrained 
system is derived easily. In Sec. lVII the quantum-classical 
stationary density matrix for a system with holonomic 
constraints is derived from the Dirac quantum-classical 
bracket. In Sec. IVIII rigorous linear response for con- 
strained systems is introduced. Finally, conclusions are 
given in Sec. IVIIII 



II. DIRAC BRACKET FOR CLASSICAL 
SYSTEMS WITH HOLONOMIC CONSTRAINTS 

Consider the phase space point X = (i?, P) , where R 
and P are coordinates and momenta, respectively, of the 
system under consideration. Let 



1 

-1 



(1) 



be the symplectic matrix, then it is well known 
that Poisson brackets can be written as 



{a,b} 



2N 



da 

2J = 1 



db 

ax" 



(2) 



where a(X) and b{X) are two arbitrary phase space func- 
tions and 2N is the dimension of phase space. The struc- 
ture of Eq. |2Jl has been exploited to introduce Hamilto- 
nian non-canonical flows [l^ . non-Hamiltonian 10. 
and constrained flows Q . Here it is summarized how to 
generalize the structure of Eq. (0) in order to define equa- 
tions of motion for systems with holonomic constraints. 

Consider a system with Hamiltonian TCq and a set of 
phase space constraints 



Following Dirac one can introduce the matrix 

C„,M&,&).^||b^.|| (4) 

t,j=l 

and its inverse (C^^)^^, where a, (3 — 1,...,2L Note 
that, when considering Poisson brackets of the con- 
straints, the convention of first evaluting the brackets 
and then imposing the constraints must be followed |^. 
By defining an antisymmetric matrix 

k,m—l a,/3— 1 

the Dirac bracket can be introduced as 

2N 



r h\ ^ da jj db 
{a,b}D= ^ jtyA 



dx, dx, 



(5) 
(6) 



Equation was originally given by Dirac in the equiv- 
alent form 

2/ 

{a, b}D = {a, 6} - E {«' iC-')^t^ b} . (7) 

a, 13=1 

Constrained phase space flows are then defined by 

d = {a^Hojo , (8) 

so that the Hamiltonian and any function of the con- 
straints is conserved because, considering an arbitrary 
function /(^a) of the constraints, it is easy to verify that 
{f{£,cr),'Ho}D = 0. The above formalism due to Dirac 
has been specialized to non-relativistic systems with holo- 
nomic constraints |^. 

In order to see how this can be achieved, consider a 
system with a number I of holonomic constraints in con- 
figuration space aa{{R}) = 0, a — l,...,l. Consider 
also the following additional constraints &a{{R,R}) = 

Y.Ziidcra/dR) ■ Pi/M = 0, a = where M are 

the particle masses. The whole set of constraints can be 
denoted as 

{£,1, ■■■,S.2l) = (ci, (7;, (71, (T/) . (9) 

Defining the matrices 



V M dR ' dRdR ' M 



i,k=l 



1 d(7a d'^a/3 P 



M dR dRdR M 



and 



Ca(X) = a = l,. 



.21 



(3) 



_ 1 dag dap 
~ ^ M dR' dR 



(10) 



(11) 
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Then one finds 



C = 



z 
z r 



and 



c-i = 



z-^rz-i -z-i 
z-i 



(12) 



(13) 



The time evolution of the phase space point under the 
constrained dynamics is then given by the Dirac bracket. 
The expUcit equations of motion are 



R = {R,no}D 



p 

M 



d(Ja 



(14) 
(15) 



where F are the forces acting on the particles. The Aq 
are the exact Lagrange multipHers 



f3=l 

' ' P d'^ap , F dap 



E 

/3=1 



M M dRdR M dR 



(16) 



III. GENERALIZED QUANTUM MECHANICS 
AND QUANTUM-CLASSICAL BRACKETS 

Consider now the realm of quantum mechanics and 
an arbitrary set of quantum variables XA, A = 1, ...,rt. 
The commutator [xa,Xn]- XaXn - XkXa {A,N = 
l,...,n) was written in symplectic form |15| using the 
matrix defined in Eq. (Q: 



[xa,Xn] = [ Xa xn]-'B' ■ 



XA 
XN 



(17) 



Considering the Hamiltonian operator of the system Hq, 
the laws of motion in the Heisenberg picture are written 
as 



dXA * r £r - 1 
— = -[i?o Xa]-B 

where the Liouville operator 



Ho 

XA 



= i^XA, 



Ho 



(18) 



(19) 



has also been introduced. 

As it was showed in Ref . [Tsf , the laws of quantum me- 
chanics can be generalized by defining the antisymmetric 
matrix operator 



C 

-C 



(20) 



with C arbitrary operator or c-number, and defining non- 
Hamiltonian quantum brackets (commutators) as 



[xa,Xn]- ^ [xa Xn] -T^- 



XA 
XN 



(21) 



Generalized equations of motion arc then defined by 

= iCxA. (22) 



dXA * r £r - 1 n 



Ho 

XA 



Quantum-classical systems are defined in terms of 
quantum variables xa in a Hilbert space which depends 
from the classical phase space point X, i.e. xa{X). Time 
evolution in the Hilbert space and in classical phase space 
are coupled consistently. The energy Eo of such systems 
is defined in terms of a quantum-classical Hamiltonian 
operator Hq = Ho{X) so that Eq = Tr' / dXHo{X). 
The dynamical evolution of a quantum-classical operator 
x{X) is given by the quantum-classical bracket which is 
defined in terms of the commutator and the symmetrized 
sum of Poisson bracket Following Refs. [lol [lll|. 

in Ref. 15] the quantum-classical bracket was easily 
casted in matrix form as a non-Hamiltonian commuta- 
tor. To this end, one can consider the operator A de- 
fined in such a way that applying its negative on any 
pair of quantum-classical operators xa(X) and xn(X) 
their Poisson bracket is obtained 



2N 



9XA „^ dx 



{xa,Xn} = -XA{X)kxN{X) = ITT^ 



N 



i,3 = i 



dx, dx. 



If one defines 



(23) 



(24) 



in Eq. (|2(J|I then a new matrix operator is introduced as 

^qc ^ ^qCQS (25) 

and the quantum-classical law of motion can be written 



dtXA 



[Ho Xa]-'D'^ 



Ho 

XA 



(Ho-.Xa) = iCxj 



(26) 

where the last equality introduces the quantum-classical 
Liouville operator in terms of the quantum-classical 
bracket. 



IV. QUANTUM-CLASSICAL DYNAMICS WITH 
HOLONOMIC CONSTRAINTS 

Consider a quantum-classical system with Hamiltonian 
Ho{X) and a set of holonomic constraints, as specified by 
Eq. (13, acting on the classical bath variables X. Then 
introduce an operator such that it can be used to 
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give the negative of the Dirac bracket of two arbitrary 
quantum-classical variables 



2N 



dxA „d9x 



N 



{xa,Xn}d = ~XAiX)A^XNiX) = J2 Qx- ' 

(27) 

where has been defined in Eq. (|3J|. Using A^, one 
can define 



D 



2i 



and the constrained matrix operator 



(28) 



(29) 



Using and , defined in Eq. ^ and ^ respec- 
tively, equations of motion for quantum-classical systems 
with holonomic constraints are defined by 



dtXA 



[Ho 



XA 



V 



D 



Ho 

XA 



{Ho,xa)d = iC XA- 
(30) 

In Equation J^OJ), the Dirac quantum-classical bracket 
{...,.. .}£) and the Dirac Liouville operator iC^ have 
been introduced. Equation (|30|l is one of the main re- 
sults of this work. It introduces the correct algebraic 
quantum-classical evolution for systems with holonomic 
constraints. In the following, all the other results of this 
paper will be derived from this equation. 

One can derive the fist consequence of Eq. H30|l by con- 
sidering the time evolution on an arbitrary function f{S,-y) 
of the holonomic constraints. This is of course given by 
the Dirac quantum-classical bracket {Ho, f{^-y))D which 
can be written explicitly as 



(ffo,/(e.))^ 



^0 7(^7 



1 

-1-^ 

1 



Ml 

2i 



2i 

Ho,f{^-f) 





Ho 







{{HoJi^y)}D 



{f{^^),Ho}D 



(31) 



Equation (|31|l shows explicitly that the Dirac quantum- 
classical bracket is defined in terms of the symplectic 
commutator (in the secon equality, the first term on the 
right hand side) of two variables minus one half of the an- 
tisymmetric combination of classical Dirac brackets. The 
symplectic commutator of the quantum-classical Hamil- 
tonian Ho with the arbitrary function 7(^7) is identically 
zero because, by hypothesis, the phase space constraints 
f^, 7 = 1, . . . , 2/, do not involve quantum degrees of free- 
dom. Then, in order to see what is the effect of the 
quantum-classical Dirac bracket on the arbitrary func- 
tion of the constraints f{C-y), there remains to be consid- 
ered the action of the classical Dirac bracket {...,.. 



For example, take 
{HoJ{^,)}d = {HoJi^^)} 

- {^o,c4(c)^a{eo,/(C7)} 
= {Ho,my)} 

- {-^0,^/3 }(C)^^C„^-^^ 



{^0,/(e7)}-E{^0'^M} 



21 



(32) 



In the same manner, one finds that {fi^-y), Ho}d = 
so that, as desired, the Dirac quantum-classical bracket 
leaves invariant any function of the holonomic con- 
straints. 



V. REPRESENTATION OF THE 
QUANTUM-CLASSICAL DIRAC BRACKET IN 
THE ADIABATIC BASIS 

Equation H3U|) is an algebraic equation. In order to ac- 
tually perform numerical calculations one needs to intro- 
duce a basis. The adiabatic basis is particularly suited 
to represent quantum-classical equations of motion 0, 
discuss quantum-classical statistical mechanics 0, ob- 
tain surface-hopping algorithms 0, 0, @ i generalize 
quantum-classical dynamics in order to have constant 
temperature dynamics on the classical bath degrees of 
freedom To define this basis, consider the following 
specific form of the quantum-classical Hamiltonian 



Ho 



K + ^R) 



2M 
+ h(R) , 



(33) 



where K is the kinetic energy operator of the quantum 
degrees of freedom and $(-/?) is the potential energy oper- 
ator coupling quantum and classical variables. It is well 
known that adiabatic states are defined by the eigenvalue 
equation 



h{R)\a;R) = Ea{R)\a;R) . 



(34) 



Before finding the representation in the adiabatic basis, 
Eq. I|30(l can be rewritten in a more explicit form 



{Ho,x)d = ^[Ho,x]-l{{Ho,x}-{x,Ho} 



/ 21 

- Yl {^o,^a}(C-i)a^{^^,x} 



- E {x,^i.}iC-')M^„Ho}\ , (35) 
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where indices with an overbar have been introduced to 
indicate constraints. The first two terms on the right 
hand side of Eq. (p?^ give the standard quantum-classical 
bracket as given if Ref. "2]. The other two terms in 
the right hand side of Eq. (|35|l pertain to the quantum- 
classical Dirac bracket and were not analyzed in previous 
works. In order to evaluate these terms, one has to re- 
mind the definition of in Eq. H13|) and the fact that 



_ fda P_ d'^CT \ 
'dR ^ laR' M ' dRdR) 



dR [ 'MdR 



(36) 
(37) 



Then consider 



I 

= ^ {7?o,a4(z-irz-i)3^{a^,x} 

a,/3=i 

-{^o,c^a}(Z"^)a^{CT^,x} 
+{Ho,&^}{Z-%-^{a-0,x} . (38) 

Recalling that {Ho^cra} = — o'g = 0, Eq. becomes 

EaJ=l {^o,ea}(C-i)a^U^,x} 
I 

= ^ {Ho,&^}{Z-^)^-0{a-^,x} ■ (39) 
a,/9=l 

Analogously, the last term in the right hand side of 
Eq. ^ is 

El,=i {x,&}(c-i)^,U.,i?o} 

= - ^ {x,f^p}(Z-^)p,R,iIo}. (40) 

fl^u—l 

Using Eqs. H38|l and (|40|l . the Dirac quantum-classical 
bracket can be rewritten as 

{Ho,x)d = ^[Ho,x]-\{{Ho,x}-{x,Ho}) 

+ 2 ^ ({^o,ap}(Z-i)p,K,x} 

fi.r>—i 

+ {x,a,}iZ-')-^,{&„Ho}) . (41) 

The last two terms in the right hand side of Eq. H41|l can 
be written more explicitly as 



1 ' 

- J2 {{Ho,a,}{Z-'),,{a,,x} 



{x,ap}(Z-i)pp{dp,iJo} 



1 ' 

- F 

2 ^ 



1 dHo 9fTp P P d^ap 



M dR dR M M dRdR 



^(y-i\ da„ die dxd(Ju ,rj-u 

^ ''^""dRdR dPdR^ '^^^ 

( P_ P_ d'^a-y 1 9gp(9go ^ 
^ I M ® M ' dRdR ~ 'M~dR~dR \ 



(42) 



With the above equations, the quantum-classical Dirac 
bracket can be written finally as 



(i?o,x)i 



xyz I 



-[Ho:X]-l{{Ho,x}-{x,Ho}) 
' 1 dHo dan P ^ P d^un 



M dR dR M M dRdR 



x(Z-i)^p 



9crp die die 9crp 



dR dR dP dR 
P P d^cTp 1 dap dHo' 



\MM dRdR M dR dR 

(43) 

Equation (|43|l is in a form suited to be represented in the 
adiabatic basis. In such a basis one can write 



{a-R\{Ho,ic)DW-R) 



/ , ^'-aa' ,I3I3'X 
130' 



where 

^iCaa',0l3'X 

I3p' 



13/3' 



1313' 



(a:R\ 



{iCaa',l3p' +lC''^Z,l3f3')X 



I3f3' 



(44) 



-[i?o,x] 



\({Ho,ii]-{x,Ho] 



\a';R) . 
(45) 

The Liouville operator iCaa'.f3/3' was given in Ref. 

i^aa',l3l3' = ii^aa' ^0/3^0' l3' + iLaa' ^a^^a' 13' — Jaa' ,1313' , 

(46) 

where UJaa' = ^~^{Ea{R) ~ Ea'{R)), 



P d 1 
MdR ^ 2 



'\ d 



dP 



(47) 



is a classical-like Liouville operator which makes 
quantum-classical variables evolve on a constant en- 
ergy surface with Helmann-Feynman forces given by 
(1/2)(F" -FF"'), and 



Jaa',f3(3' = —5a'l3'da(3 
- 5al3d*^,p, 



p h 



p h 



d_ 
Wl 
d_ 
'dP 



(48) 
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is an ofF-diagonal operator, realizing nonadiabatic tran- 
sitions, which is defined in terms of the nondiabatic cou- 
pling vector dap — {a\d/dR\P). The operator iC^^a' ,i3p' 
which imposes the constraints in the quantum-classical 
dynamics is 

E-rcon 13/3' 
(30' ^'-aa' ,I3I3'X 

1 V-/ 1 dHodn-^ P _ P 9V„ ' 



M dR OR M M dRdR 



( P P d'^ap 1 ScTp dHa 



M M dRdR M OR OR 



\a-R) 



(49) 



By defining F°'P = -{a: R\dHo/dR\l3; R) and F" = 
—dEa{R)dR, using F"'^ = F^Sap + fiujapdap, adding 
and subtracting the term 



1 , ' 



1 dUr, 



^^(Z-^)p.^^<5„^<5„,,, (50) 



and doing some algebra one finally obtains 



■ /-con 



dR dP 



I30a'f3' 



1 

~2M ^ {TT'^al3dap5a'j3' 

fi,L> — l 

+ hUa'l3'da,p,5ap) ■ "^(Z"^)/iP 



d(Ta,„_^ dap d 



dR dP ' 
(51) 



where the A"" , which are the quantum-classical La- 
grange multiphers on the energy surface (l/2)(£'a + i?Q'), 
are defined as 



— 1 



M M dRdR 



The first term on the right hand side of Eq. (|51|l defines 
a diagonal operator 



{\/2){Ea -\- Ea'). The other two terms in the right hand 
side of Eq. (|51|l defines an off-diagonal operator 



Tcon 

'Jaa',PP' 



1 ' 

— ^ (htOapdapSa' p' 



2M 



huJa'p'da,p,Sap) 



da 



das d 



dR ^ '^^ dR ■ dP 
(54) 



which couples the constrained dynamics to quantum 
transitions between the adiabatic states. Its effect is, on 
the one hand, to modify the probability that quantum 
transitions take place and, on the other one, to realize 
the back-reaction of the quantum transitions on the con- 
strained momenta. It is convenient to finally cast the 
quantum-classical Dirac-Liouville operator of Eq. H44|) in 
the following form 



iC 



D 

aa',PP' 



[iuJaa' + iLaa') ^apSa'P' " J, 



D 

aa' ,PP' 1 



(55) 



where 



P d 

MdR ' 2 



_d_ 

dP 



' dac d 
'dRdP 



(56) 



and 



jD 7 1 rcon 

'^aa',PP' — •Jaa',Pp' + •Jaa',PP' 



P 

M 



lap 



2 (£■ • dap) 



M "'"''3' 



^+2 



ft i^a'P'd*a>p' 



-d* 



d_ 

dP 

d_ 
' dP 



where the constrained frequency is 
n I . . ^ ^ da 



•^ap 



^ap 



dajy 



^ M dR^ ^^"^ dR 
/i,i/=i / 



Sa'p' 
SaP , 

(57) 
(58) 



The time evolution of any dynamical variable is given 
obviously by 



(59) 



PP' 



aa'dai, d 



dR dP ' 



(53) 



A. MOMENTUM-JUMP APPROXIMATION 
FOR CONSTRAINED SYSTEMS 



whose action is to enforce the constraints while time evo- 
lution takes place adiabatically on the energy surface 



When performing numerical calculations on systems 
with many degrees of freedom 0, |^ IT^ [l7| the action of 
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the operator in Eq. H48() is usually evaluated within the 
momentum-jump approximation ^18j . This approxima- 
tion can be derived easily for the operator in Eq. H57|l 
defined in terms of the constrained frequency ujPa in 
Eq. (|58|l . In analogy with what shown in Ref. [igj in 
order to derive the momentum-jump approximation for 
constrained system, one can first consider one of the two 
terms in the right hand side of Eq. (|57|l 



from that given in Eq. H64() is obtained. Instead, Equa- 
tion (|64|l . with the momentum-jump operator for the con- 
strained system in Eq. H()2I) . was derived correctly from 
the fundamental Dirac quantum-classical formalism and 
it provides the correct formula for constrained systems. 



VI. STATIONARY DIRAC DENSITY MATRIX 



h f P \ d 



M 



d 



dP 



d{P ■ do/3)2 



(60) 



where dap denotes the unit vector along the direction of 
dap- Then perform the approximation 



1 + huj^i^M 



d 



d{P ■ daf3) 



exp 



d 



d{P ■ d„/3)2 



(61) 



Write the quantum-classical Dirac density matrix of a 
constrained system as pb{X). The average of any oper- 
ator X can be calculated from 

(X) = Tr' j dX pDxit) - Tr' J dX pn exp (iC^'t) x ■ 

(65) 

The action of exp (iC^t) can be transferred from x to 
Pd by using the cyclic invariance of the trace and inte- 
grating by parts the term coming from the classical Dirac 
brackets. One can write 



and write the operator J^^, ^p, as the operator J^^, 
evaluated in the momentum-jump approximation 



jD,M-.J 



P 

M 



dai3 exp 



T7 • ) exp 



d 



d{P-daf,)\ 
d 



The action of any of the operators on the right hand 
side of Eq. (|62|) on an arbitrary function of momenta is a 
translation (or jump) of the momenta itself. For example. 



exp 



d 



d[P ■ dapY 



f[P) = /(P + AP) 



(63) 



with 



AP = sign{P-dap)^{P-da0Y + nMuj^p-{P-dafi) 



(64) 



Equation H64|l illustrates the effect of the momentum- 
jump operator on a function of constrained momenta. It 
differs from that given in Ref. 'iS"] and used, for example, 
in Refs. 4, 5, 16J because of the constrained frequency 
uo^p defined in Eq. (|55|l . This result, which should have 
been used in the nonadiabatic calculations of Ref. 
is reasonable because constrained momenta cannot be 
scaled as unconstrained momenta. As a matter of fact, if 
one follows the empirical procedure of applying the un- 
constrained momentum-jump of Ref. [T^ . as it was done 
in Ref. , and uses the RATTLE procedure to im- 
pose the constraints on momenta then a result different 



Ho, ■ 



1 



{Hf), ■ ■ ■}d — {•••, Ho} 



D 



(66) 

In this equation the classical Dirac bracket terms are 
written 



2N 



Saf3 ■ 

(62) 



{Ho, ■ ■ - jo - {■ ■ ■ ,Ho}d — 

d.. 



dHo 

dX, dX, 



odHo 



dX, dX 



(67) 



When integrating by parts the right hand side, one ob- 
tains terms proportional to the compressibility 



2N 

E 



dB^dHo 
dX, dX. 



(68) 



Using Eqs. Q and ifT^ . in analogy with Ref. one 
finds easily that 



D 



— In detZ 

dt 



(69) 



Because of the compressibility in Eq. 1)69(1 . it turns out 
that the Dirac quantum-classical Liouville operator is not 
hermitian 

(iC'^y ^ -tC"" - . (70) 

The average value can then be written as 

ix) = Tv' J dX xexp [-(i/:^ + K{^)t] PD . (71) 

The quantum-classical Dirac density matrix evolves un- 
der the equation 



d_ 

dt 



PD 



Ho, Pd 



PD 



(^{Ho, Pd]d - {pd,Ho}d 

(72) 
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The stationary density matrix p^e is defined by 



iC^PDe + Pd 







(73) 



To find the exphcit expression of p/^e one can follow 
Ref. , expand the density matrix in powers of h 



PDe 



n=0 



(74) 



and look for an explicit solution in the adiabatic basis. 
In such a basis the Dirac-Liouville operator is expressed 
by Eq. (|55|) and the Hamiltonian is given by 



2M 



(75) 



Thus, one obtains an infinite set of equations correspond- 
ing to the various powers of fi 







■ T D {n)Q.a n (n)oLOi 



As shown in Ref. Q , in order to ensure that a solution can 
be found by recursion, one must discuss the solution of 
Eq. H77|l when calculating the diagonal elements p^^"" 
in terms of the off-diagonal ones p^g"" • To this end, 

USmg p^g — l^p^^ ) , Jaa,0l3' " '^aa,P'P '^'^^ 

fact that J^Q, = when a real basis is chosen, it is 
useful to re-write Eq. (|77|l in the form 



(*Lf„ + Op 



/3>/3' 



De 



2TZ 



( jD in)f}p' 
Yaa,l3f3'PDe 



(78) 



One has ^ {-iL^a " ^^o*)^ = «^aa- The right hand side 
of this equation can expressed by means of the classical 
Dirac bracket in Eq. ©. It follows that Hl^ and any 
general function /(Hq) are constants of motion under 
the action of iL^a- Because of the presence of a non-zero 
phase space compressibility, integrals over phase space 
must be taken using the invariant measure |a, l20l | 



dM = exp{-WD)dRdP , 



(79) 



where wd = J dtn^ = defL is the indefinite integral of 
the compressibility. To insure that a solution to Eq. H78|l 
exists one must invoke the theorem of Frcdholm alter- 
native, requiring that the right-hand side of Eq. H78|) be 
orthogonal to the null space of {iL^aV US- The null- 
space of this operator consists of functions of the form 
/{Hq), where /{Hq) can be any function of the adiabatic 
Hamiltonian Hq. Thus the condition to be satisfied is 



dM J2 2^ 



( jD Mf3l3' 



) f{H^) = . (80) 



To this end, there is no major difference with the proof 

given in Ref. 0: 27^ ( Jf„,^0,p^f and f{H^) are re- 

spectively an odd and an even function of P; this guar- 
antees the validity of Eq. H8U|) . Thus, one can write the 
formal solution of Eq. H78(l as 

PDe = (*-^aa + '^O ) 2^^^ [J aa^Pfi' Poe 

and the formal solution of Eq. H77|) for a ^ a' as 



(81) 



POe 



D 



nx («)aa' 
'^0 )PDe 



En 



(n)l3l3' 



(82) 



PI3' 



Equations H81() and H82|l allows one to calculate p^" to 
all orders in h once p^g"" is given. This order zero term 

is obtained by the solution of {iL^a~^ ''^o)PdI°^°' ~ 0. All 
(76) higher order terms are obtained by the action of Eaa' , the 
imaginary unit i and J^a'p/S' (involving factors of c?qq', 
^aa'' ^iid derivatives with respect to P. 

One can find a stationary solution to order h by consid- 
ering the first two equations of the set given by Eqs. H7t)|) 
and lf77|l : 



J^O,PDe 

Ho,P^nl 







(n = 0) 



-^[{Ho,P^^I}d-{p'ZHo}d 



= (0) 



2 I'^O ; PDe\ + 



(n = 1) . 

For the 0{h^) term one can make the ansatz 



POe 



(83) 



(84) 



(85) 



where Q is 



Q ^ J2 fdM 6(^)6 {C-H^) (86) 

and S{^) is a compact notation for 
I I 

= n ^('^'^) n ^('^a) • (87) 

The following expression for the order h term is obtained 



.(l)a/3 
POe 



(0)/3 



f3 



Ep — Ea 



(88) 



for the 0{h) term. 

Equations H85|l and 1)88(1 give the explicit form of 
the stationary solution of the Dirac Liouville equation 
up to order 0{Ti). This stationary solution must be 
used when calculating equilibrium averages in quantum- 
classical systems with holonomic constraints on the clas- 
sical variables. 
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VII. LINEAR RESPONSE THEORY 



then iCfit) = -T{t)iC^ and Eq. ^ becomes 



Consider a perturbed quantum-classical Hamiltonian 
H{t) =Ho + Hi{t) =Ho- AJ^{t) . (89) 
Define the perturbed Liouville operator 



{Hi{t),...}D~{...,Hi{t)}D] .(90) 



In general the perturbation can bring an additional term 
to the compressibility of phase space 



Using Eqs. lO and ((T^ . one finds 



'dlndetZ J-^ dZ^f 



dR 



l/,fl—l 



dR 



OA 



(92) 



In this case, the pertubed Liouville operator is not her- 
mitian 



where 



(z£f(t))t = -^Cf{t)-^[kf{t),...l 
...,...]+ denotes the anticommutator 



(93) 



[fi?,PD]+ = [fif Pd] 





1 






PD ] ■ 


1 




PD 



af PD 



'PDkf 



(94) 



The constrained evolution of the density matrix is ob- 
tained from Eq. H72|) by replacing /fp by H{t) 



d_ 

di 



PD{t) 



H{t),pD{t) 



+ ^[{H{t),pDmD-{pD{t),Ho}D 



(95) 



Assuming that the perturbed density matrix is poit) = 
Pdc + ^PDit) and that the system was in equilibrium in 
the distant past, linear response theory gives 



ApD(t) = - / dT&M~lC''\t^T)]iC'^\T)pDe 
J — oo 



(96) 



If one defines 



l({i(t),...}^ -{..., ijz.) , (97) 



Api)(t) 



dTT{T) exp[-i£^t(t „ T)]iC'^^pD. 



A yDe 



PDe 



^[KA.PDe\ + 



The non-equilibrium average of any quantum-classical 
operator B{X) can be calculated over the density matrix 
poit), B = Tr' J dXB{X)pD{t) in order to determine 
the response of the system to the external force. 

AB{t) = Tr' j dXB{X)ApDit) 

= f dT$B^(i-T)^(r) , (99) 
<y —oo 

where the response function is defined as 



^BA{t) 



Tr' / dXBiX;t)[iC^pDe + 7;{ 



■^i'^AiPDe 



(100) 



Equation (|100|l gives the response function for quantum- 
classical systems with holonomic constraints on the clas- 
sical variables. In analogy with the purely classical case, 
analized in Ref. (^j, if A{X) depends from the momenta 
P the response function contains two contributions in 
addition to the expression of $B^(i) given in Ref. 0. 
One of this contributions arises evidently from and 
the other comes from the factor detZ contained in the 
expression for poe- 

Because a perturbation operator which depends only 
on particle positions is usually adopted in order to derive 
expressions for quantum-classical rate constants (typi- 
cally the Heaviside function), the rate formulas derived 
in Refs |Q H, El also applies to a constrained system with 
the exception that the correct constrained stationary ma- 
trix, as derived in Sec IVII must be used. However, dif- 
ferent calculations and different perturbation operators 
may require the evaluation of the additional terms of the 
response function here derived and one must be aware of 
their existence. 



VIII. CONCLUSIONS 

In this paper the theory of quantum-classical systems 
has been generalized in order to treat rigorously situa- 
tions where the classical degrees of freedom must obey 
holonomic constraints. The formalism here presented 
has been obtained by unifying the classical Dirac bracket 
with the quantum-classical bracket in matrix form. In 
this way, a Dirac quantum-classical formalism, which 
conserves the constraints exactly, has been introduced 
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and then used consistently to formulate the dynamics 
and the statistical mechanics of quantum-classical sys- 
tems with holonomic constraints. A first result has been 
the derivation of the correct momentum-jump approxi- 
mation which takes into account that, when a quantum 
transition occurs, the momenta cannot cannot be scaled 
as if they were unconstrained because, instead, the holo- 
nomic constraints must be satisfied. Moreover, the Dirac 
quantum-classical bracket allows one to derive easily lin- 
ear response theory. This has shown that the rigorous 
response function of constrained systems contains non- 
trivial terms, which were already noted by this author in 
classical mechanics, arising from the action of the per- 
turbation operator on the phase space measure of unper- 
turbed constrained systems and from the compression of 
phase space which may be caused by the perturbation it- 
self. These terms are zero if the external perturbation is 



coupled only to the position coordinates of the classical 
degrees of freedom. 

If one considers with a wider perspective this work and 
that of Ref. TS], dealing with the introduction of non- 
Hamiltonian commutators in quantum mechanics, it can 
be realized that a unified formalism for defining general- 
ized dynamics in quantum-classical systems is now avail- 
able. There are reasonable expectations of employing in 
the future specific forms of such generalized dynamics in 
order to attack the problem of long time numerical inte- 
gration of quantum-classical dynamics. 
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